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The dynamics of a decohering two-level system driven by a suitable control Hamiltonian is studied. 
The control procedure is implemented as a sequence of radiofrequency pulses that repetitively flip the 
state of the system, a technique that can be termed quantum "bang-bang" control after its classical 
analog. Decoherence introduced by the system's interaction with a quantum environment is shown 
to be washed out completely in the limit of continuous flipping and greatly suppressed provided 
the interval between the pulses is made comparable to the correlation time of the environment. 
The model suggests a strategy to fight against decoherence that complements existing quantum 
error-correction techniques. 
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I. INTRODUCTION 



The design of strategies able to protect the evolution of a quantum system against the irreversible corruption 
due to environmental noise represents a challenging conceptual issue. In particular, since maintaining quantum 
coherence is a crucial requirement for exploiting the novel possibilities opened up by quantum parallelism, practical 
implementations of quantum computation and communication proposals require methods to effectively resist the 
action of quantum decoherence and dissipation Roughly speaking, two classes of procedures are available to 

00 ■ overcome the decoherence problem: either passive stabilization or active manipulation of the quantum state. The 
0^ ' first kind of solutions, recently formalized as Error- Avoiding Codes Q, relics on the existence of a subspace of states 
that, owing to special symmetry properties, are dynamically decoupled from the environment. The second approach, 
pionereed in j^] and closer in spirit to quantum control theory |j]|7| , embraces today a variety of sophisticated schemes 
known as Error-Correcting Codes |q- |lq] . Basically, loss of information is corrected by monitoring the system and 
conditionally carrying on suitable feedback operations. 

In this work we investigate a third strategy for reducing noise and decoherence. This strategy, which can be 
termed quantum "bang-bang" control after its classical analog ]l9| ], works by averaging out the unwanted effects of 
the environmental interaction through the application of suitable open-loop control techniques on the system. The 
basic idea is that open-system properties, specifically decoherence, may be modified if a time-varying control field 
acts on the dynamics of the system over time scales that are comparable to the memory time of the environment. In 
particular, we work out an exact model for a two-state quantum system (qubit) coupled to a thermal bath of harmonic 
oscillators, where decoherence is dynamically suppressed through repeated effective time-reversal operations on the 
combined system + bath. Although the phenomenon is mathematically reminiscent of the quantum Zeno effect [ 20 1 , 
the essential physical idea comes from refocusing techniques in nuclear magnetic resonance spectroscopy (NMR) [21|. 



Since the discovery of spin echoes in 1950 22 1, clever pulse methods have been developed in NMR to eliminate much 
of the dephasing arising from variations in the local magnetic field acting on each spin. Because the latter effect can be 
thought in terms of an interaction with some classical environment, it is not obvious a priori whether similar techniques 
work in the presence of a quantum mechanical environment and purely nonclassical effects like entanglement. Our 
result answers this question in the affirmative, and points out the role of the reservoir correlation time as a further 
parameter to be engineered in the struggle for preserving quantum coherence. 

The plan of the paper is the following. In Section II, a general model of a two-state system interacting with a 
thermal environment is reviewed and its decoherence properties in the absence of control recalled. In Section III, the 
evolution under the action of a sequence of perturbative kicks is analyzed. Complete quenching of decoherence is 
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established as a limiting situation. In Section IV, the conditions for an effective decoherence reduction are clarified 
with reference to a variety of possible environmental configurations and the method is compared to different quantum 
error-correction techniques. We close by discussing possibilities for future work. 



II. SINGLE-QUBIT DEPHASING MECHANISM 



Our goal is to investigate how decoherence properties of an open quantum system may be modified through the 
application of an external controllable interaction. Decoherence is a process whereby quantum systems lose their ability 
to exhibit coherent behavior such as interference |p3|-p6f . We start by introducing a model that allows investigation 
of the problem in its simplest nontrivial configuration. The physical system we are interested in is a single two-state 
quantum system, representing the elementary memory cell of quantum information (qubit). Although not strictly 
necessary, it will be convenient to think of the physical qubit as realized by a spin-I/2 system, which will provide us 
with direct reference to the well established language of Nuclear Magnetic Resonance pM and the rapidly growing 



field of NMR quantum computation []27 28 1. Decoherence arises due to the coupling to a quantized environment, here 
schematized as a continuum of harmonic modes. We assume that the dynamics of the overall qubit + bath is ruled 
by the following Hamiltonian: 

H = H S + H B + H SB = huj ^+J2 ^kb{b k + Ko z (g k b{ + g* k b k ) , (I) 

k k 

where the first and second contribution Hg and Hg describe, respectively, the free evolution of the qubit and the 
environment, and the third term Hsb describes a bilinear interaction between the two. a z is the standard diagonal 
Pauli matrix, with qubit basis states denoted as i = 0, 1, while b\, b k are bosonic operators for the fc th field mode, 
characterized by a generally complex coupling parameter g k . In the Schrodinger picture, the state of the combined 
system (S + B) is represented by a density operator ptot(t) and the reduced qubit dynamics is thereupon recovered 
from a partial trace over the environment degrees of freedom: 

Ps{t)= E Pij(t)\^ij\ = ^B{Ptot(t)} ■ (2) 
i,j=0,l 

Hamiltonian ([!]), which corresponds to a special case of the so-called spin-boson problem has been used by many 
authors to model decoherence in quantum computers ^j-Q. In particular, we adhere closely to the notations of [gj. 
The basic fact about the dynamics induced by ([l]) is that, since [a z , Hq] — 0, the interaction with the environment has 
the two memory states |0), |1) as eigenstates. In other words, the model describes a purely decohering mechanism, 
where no energy exchange between qubit and bath is present. In NMR terminology, this implies that no Ti-type of 
decay takes place plj. Equivalently, in terms of errors, only phase errors are introduced. However, neglecting the 
effects associated to quantum dissipation is justified, in a sense, by two related reasons: energy exchange processes 
not only produces amplitude errors which need to be corrected even in the classical computation, but they typically 
involve time scales much longer than decoherence mechanisms. In addition, being exactly soluble, the model (|) has 
the advantage of allowing a clear picture of the decoherence properties in the absence of control. To this end, since 
spin populations are not affected by the environment, the relevant quantity is the qubit coherence poi(t) (of course, 
Pio(t)=P* i(t))- 

It will be convenient to move to the interaction picture associated to the free dynamics (H$ + Hb), corresponding 
to the transformed state vector 

p tot (t) = e Ws+H B )t/n ptot{t) e - i{Hs+HB)t /n (3) 

and to the effective Hamiltonian 

H(t) = Hsb (i) = ha, £ (^e^* + ^e"^') . (4) 

k 

Time evolution is determined by the time-ordered unitary operator 

Utot{t Q ,t) =Texp| - ijT dsH(s)^ , (5) 
which can be evaluated exactly and can be written, up to a global c- number phase factor, in the following form: 
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^ Ot (to ) t) = exp|^^(6^*°4(*-*o)-^e-^ to ^(t-to))| , (6) 



k 

where the following complex function has been introduced: 

Z k (M) = ^(l-e i »* At ) . (7) 
uj k V ) 

Note that the separate dependence of (Q) on both initial time t and evolution interval (t — to) is consistent with 
the time-reversal property Utot(t,to)Utot(to,t) = 1. A nice discussion on the entanglement generated by Utot(to,t) 
between qubit and field states is given in jj| . We are interested in calculating 

An(t) - (0\TT B {U tot (t ,t)p tot (t )ul t (to,t)}\l) . (8) 

This can be done without approximations under two standard assumptions: 

(i) qubit and environment are initially uncorrelated, i.e. 

Ptot(to) = As (*o) ® /Ob (to) ; (9) 

(ii) the environment is initially in thermal equilibrium at temperature T, i.e. 

pB{h) = J] PB, fe (T) = [] (l - e"*"*) e-^ b W . (10) 

fc fc 

In Eq. (|Io|), /3 = l/feT. being fc^ the Boltzmann constant. For simplicity, we choose henceforth units such that 
ft = ks = 1. Conditions (|)-(10) are easily translated to interaction picture and, inserting (|^) into (|8|), the problem 
is reduced to a single-mode trace: 

Poi(t) - Poi(to) • 1] Tr 4 PbAT) V{^ ta i k {t - t )) } = An (to) • e^* '* , (11) 

k 

where, in the first equality, the harmonic displacement operator X>(£fc) is given by 

2>(&)=e*te- 6 »«, (12) 

and the second equality defines the time-dependent function ro(to,t). The final step is to recognize that, for each 
mode, the quantity in curly brackets in ([ll]) is nothing but the symmetric order generating function for a thermal 
harmonic oscillator Thus, the explicit expression for ro(to,t) is 

r (t , t) = T (t - to) = E la V° )|2 -th . (13) 

Being To a real function, it correponds to pure damping in (pd|). Accordingly, this function characterizes completely 
the dynamics of the decoherence process destroying the qubit phase information. Of course, the complete evolution 
in the original Schrodinger picture includes oscillation at the natural frequency ujq, 

Aoi(t) = e iw °( t - t °)- r °( t - t °) aoi (to) , (14) 

while, using (||), one can check that pa(t) = pu(to), i = 0,1. 

Deeper insight into the time dependence of the decoherence process is gained if the continuum limit is made explicit 
in (13). By substituting |£fc(t — to)| 2 through (Q), we get 



/■oo r 

r (t-to)=/ du y^5(u - u k )\g k \ 
Jo L k 



. Ul \ 1 -cosw(t-to) 

4coth f2r) ^ ■ (15) 



The quantity in square brackets is known as the spectral density I(u>) of the bath. It turns out that, once the initial 
state is specified, complete information about the effect of the environment is encapsulated in this single function. As 
a general feature, the spectral density is characterized by a certain ultraviolet cut-off frequency to c such that I{u>) — > 
for uj > u) c . Although the specific value of uj c depends on a natural cut-off frequency varying from system to system, 
the existence of a finite ui c is always demanded on physical grounds. Indeed, assuming that the environment does not 
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have a high frequency cut-off means, generally, that energy can be dissipated instantaneously (3^]. If, for instance, 
decoherence arises from the coupling to a phonon field, the natural cut-off frequency can be identified with the Debye 
frequency. In general, r c ~ ui^ 1 sets the fastest time scale (or the memory time) of the environment. When a specific 
choice will be useful, we will assume a spectral density with the following functional form: 

J(w) = ^"e-^, (16) 

the parameter a > measuring (in suitable units) the strength of the system-bath coupling and the index n > 
classifying different environment behaviors P, p9|]3l| ] ■ 

In addition to r c , another time scale rg ~ T^, associated with the temperature of the bath, is expected to play a 
major role in the evolution of the qubit coherence. This is manifest if Eq. ( [l5| ) is written in the equivalent form 

r (t - to) - 4 J™ dcj 1(c) (2fS(u, T) + l) l-cos^t-M 

where n(tu,T) — exp(— w/2T)cosech(w/2T) is the average number of field excitations at temperature T. In this way, 
effects due to the thermal noise are formally separated from the ones due to purely vacuum fluctuations. Because 
of the different frequency composition, the two kinds of fluctuations dominate on different time scales, the relative 
importance of vacuum to thermal contributions being determined by Tp. This multiplicity of time scales is one of the 
factors that makes the decoherence process quite complicated. It should be important to keep in mind that there 
is no generic decohering behavior, and in particular that the qubit dynamics depends crucially on both temperature 
and the details of the spectral function I(u>). We will comment further on this point later. 



III. PULSED EVOLUTION OF QUANTUM COHERENCE 



Let us now define a procedure aimed at modifying the decoherence properties discussed so far. We choose to 
implement it as a suitable perturbation acting on some observables {Oi} of system S. This is obtained by adding to 
(|l|) a time-dependent term ji(t)Oi, where the input functions {7«(i)} are assumed to be programmable at will, 
e.g. a given schedule of time- varying magnetic fields in the case of a spin-qubit. In control terminology, this realizes a 
so-called open-loop configuration |jj . Closed-loop (or feedback) configurations have been also proposed in recent years 
to manipulate decoherence in some quantum optical systems [ ^2[ . There are many possible choices for the control 
Hamiltonian. Leaving aside the problem of controllability of the system (|l|) at the abstract level, we make here a 
pragmatic choice, partly suggested by semiclassical considerations. 

Suppose that the perturbation we add is able to induce spin-flip transitions. By inspection of the spin-bath 
interaction Hamiltonian Hsb, opposite contributions arise when the spin belongs to the down or up eigenstate. Since 
a relative minus sign will be present during time evolution, the effect of the H$b coupling will eventually average out 
provided the spin is flipped rapidly enough. There are two reasons why a mechanism of this kind is expected to be 
possible. First, as we already mentioned, similar methods are routinely used in NMR experiments to get rid of the 
effects of some unwanted interactions For instance, in the so-called spin-flip narrowing, line broadening resulting 
from magnetic dipolar coupling is reduced by repetitively flipping the spins among distinct energy configurations. 
The major difference here is that the undesired decohering coupling is contributed by the infinitely-many quantized 
degrees of freedom of the heat bath. The second reason is related to the finite response time r c of the environment 
itself. In general, if perturbations act on the system more rapidly than this fastest time scale present in the bath, 
we expect that memory effects would be relevant, and that this non-Markovian dynamics would eventually lead to 
modified decoherence properties. 

Having anticipated the intuitive idea, we start by writing the total Hamiltonian as 

H(t) =H +H rf (u ,t) , (18) 

where H r f(u>o,t) represents a monochromatic alternating magnetic field applied at resonance. At this level, the 
problem can be related to variants of the complete spin-boson model, none of them is exactly solvable f29|| . The 
main simplification we introduce is to substitute the continuous-mode operation, where the actions of the bath and 
the controlling field are necessarily simultaneous, with a pulsed-mode operation analogous to classical "bang-bang" 
control. Then, under the working assumption that the typical decoherence time of the system and the duration of 
pulses define two widely different time scales, we solve the problem by putting 

(i) Hsb = within each pulse; 

(ii) H r f = between successive pulses. 
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As usual, we invoke the rotating wave approximation and only look at the co-rotating component of the rf-field. The 
radiofrequency perturbation is assumed of the following form: 



with tp^ = to 



H rf (uo,t) = J2 V {n \t) [cos L (t - t ( p l) )y x + sin foj (t - 4 n) ))< 

n=l 

nAt, n = 1, 



(19) 



, np , and 

V {n) {t) = 



V 




4 n) < t < 4 n) - 

elsewhere 



tp , 



(20) 



Eqs. (|T^)-(20) schematize a sequence of np identical pulses, each of duration Tp, applied at instants t — tp. The 
separation At between pulses is assumed to be an input of the model. The amplitude of the field is equal to V during 
each pulse and will be further specified below together with Tp . In order to depict the evolution associated to a given 
pulse sequence, it is convenient to think the latter as formed by successive elementary cycles of spin-flips, a complete 
cycle being able to return the spin back to the starting configuration. We begin by analyzing the time evolution 
during the first spin cycle. 

A. The elementary spin-flip cycle 

As in Sec. II, we exploit the interaction representation (^|). Some instructive steps of an alternative derivation 
based on Heisenberg formalism are sketched in the Appendix. The interaction picture transformed Hamiltonian is 
now 

H(t) = H SB {t)+H rf {t) , (21) 
where Hss(i) is given in (Q), and H r f(t) is evaluated by using ( |ig| ) and the properties of Pauli matrices. We get 



H rf (t) = V V {n) (t) e M0 *'p > a x e- luJ0 ^ t{ r 



(22) 



n=l 



According to (p2|), the time dependence due to the rotating field is completely removed within each pulse. In fact, the 
interaction representation (||) on spin variables is identical, at resonance, with the description in the rotating frame 
associated to ([l9]) [^l). The counter-rotating term that is omitted within RWA is seen to be negligible at resonance. 
For the first spin cycle, np — 2 and we have the following sequence: evolution under HsB{t) during t < t < tp* 1 ; 

pulse Pi at time tp^; evolution under HsB(t) during tp + Tp < t < tp] pulse P 2 at time tp\ After a total time 
ti = to + 2A< + 2rp, the first cycle is complete. In terms of evolution operators, we have: 



Utot(to, t P \ 



U P (t ,tx) = Up 2 U Pi [Up 1 1 Utot{t i p ) +T P ,t i p)U Pl 
We can read the evolutions in the absence of rf-field directly from (^), for instance 



U tot (^+rp,t^)^e Xp \^Y,( b l e 



a iu> fc (to+At+Tp) 



Z k (At) -b k c 



-iui k (t +At+Tp) f* 



(23) 



(24) 



Concerning the evolution operator associated to the generic j^ 1 pulse, this is found by exponentiating (|22|) (n = j): 



U P] =exp j -iVrpe^^p a^-^^p \ 

In order to proceed, we have to say more about pulses. We require 

{U Pl ,(7 z } = 0, [Up 2 U Pl ,a z ] =0, 



(25) 



(26) 



{ , } and [ , ] denoting anticommutator and commutator respectively. These conditions imply, as expected from the 
intuitive explanation, that P\,P-2 are 7r-pulses, satisfying IV Tp = ±7r. To simplify things, we imagine that H r f(t) is 
large enough to produce (almost) instantaneous spin-flips. Accordingly Tp — > henceforth and we have to deal with 
ideal "kicks" of infinite power. Putting things together, we find 
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Up 2 U Pl = _ e ^o^(ti-to) ^ (27) 



U^U tot {tf ,tf)U Pl = expj - fe^+^At) - bke -^o+At) CkiAt fj 

^ k 

and we are therefore ready to write down the final result for the cycle evolution operator: 



(28) 



U P (to,h) = expj^ ^(h -to) + y ^ (bU luJkt "Vk(At) - & fee -^*°^(At)) j , (29) 



k 

where, as before, c- number phase factors have been omitted and 



% (Af) = £ fc (Ai)(c^ At -l) = _^i(i_e^ At ) 2 . (30) 



It is interesting to compare the evolution described by (|2^) with the one in the absence of pulses. By recalling (||), 
evaluated at time f i = to + 2At, we report two differencies: a phase factor proportional to a z and the duration of the 
cycle; a combination rjk(At) oc (£k{At)) 2 in place of £fc(2Ai). The first difference corresponds to the fact that, due 
to the pulses, the oscillation at the natural frequency u>o is lost once the evolution is transferred back to Schrodinger 
picture (see (|14])). The second difference, as it will be seen in a moment, is the signal that decoherence properties are 
modified. 



B. Decoherence properties after a pulse sequence 

p is to general] 

ending at time 



The next step is to generalize the description to an arbitrary number N of elementary spin-flip cycles, the rfi 1 cycle 



t n =t + 2nAt, n = l,...,N, (31) 

and the number of 7r-pulses involved in the sequence being rip = 2N. This is straight forward since Eq. (|2^) enables 
us to write down the evolution operator for the nfa cycle: 

Up(tn-i,t n ) = expjiwo y2Ai + (ble^-^kiM) - ^^-^(Ai)) j . (32) 

The time development corresponding to N cycles is then governed by the time-ordered finite product 

U { p N) (to,...,t N ) = Up(t N -i,t N )...Up(h,t 2 )Up(to,t 1 ) . (33) 

Note that, at variance with conventional spin-flip narrowing U p N ^ ^ {Up) N , the dependence on intermediate 

times {tj} in the sequence being introduced by the environment dynamics. A closed formula for U P \{tj}) can still 
be found quickly since, by neglecting state-independent global phase factors that are irrelevant to density matrix 
propagation, we are allowed to treat the factors in (R3|) as commuting operators. We get 



where 



Up N \t Q , At) = exp y 2iVAi + y ^ (b\e^ to m (N , At) - he^^UN, At)) J , (34) 



N 

Vk (N, At) = Vk (At) e 2i( - n -^ At , (35) 

n=l 

and definition ( |3l|) has been exploited. Of course, the results of Sec. HIA are recovered for N = 1. A more interesting 
check can be done by relating the evolution (^54|) to the corresponding propagator in the absence of pulses. The trick 
is to switch back the minus sign in the definition (|30| ) of r)k(At), that has been recognized as the key dynamical effect 
due to the pulsing procedure. Since ^(Ai)(e 41i ' fcA M^ 1) = £fc(2At), we are left with 



G 



N 



&(2At) 2 e 2 ^"- 1 ^^ 4 = £ fe (2iVAt) = e*(tjv - *o) (36) 



71=1 

in place of (|3^). Therefore, this procedure gives 



U tot (t , t N ) = cxp |^ ( b l^ kto tk(t N - t ) - b k e-^C k (tN - to)) | , 



(37) 



which does agree with the direct evaluation based on (|^). It may be worth some times to let the connection with 
unperturbed evolution be explicit, which is done by writing 



N 



r lk (N,At)=Z k (2NAt)-2Z k (At)J2e Mn ~ 1) " kAt , (38) 

n=l 

where the contribution due to the pulse sequence shows up in the form of a typical interference factor. 

The decoherence properties corresponding to the pulsed qubit + bath evolution ( 34 ) are derived through the steps 
outlined in Sec. II for the unperturbed case. Incidentally, we still expect that spin populations are unchanged after a 
sequence of N complete spin cycles, while qubit coherence at final time t^ — to + 2N At is given by 



The result is 



where 



Poi{t N ) = (0\Tr B {U { p N Hto, At)ptot(to)U ( p N) Hto, &t)} \1) . (39) 
Pm(t N ) = e--o(^-to)-rp(JV,A 4) > (40) 

r^j^JaM^te). (4i) 

Comparing with ([l3|), the mathematical prescription for decoherence in the presence of pulses looks very simple: just 
use rj k (N, At) instead of £fe(t/v — 1 ) for each mode of the bath. However, the final effect leading to ( |4l|) is not easy to 
figure out and it is useful to compare first the decoherence due to a single mode with frequency ui with and without 
pulses respectively. Apart from identical time- independent factors, we have to consider 

N-l 

\ri(N,ujAt)\ 2 =4(l-coswAi) 2 N+ ^ 2ncos (2(N - n)u)At) , (42) 

n=0 

vs. 



\£(N,ujAt)\ 2 = 2(1 - cos 2NwAt) = 2(1 — cosw(£jv — to)) , (43) 

where (0) and (|35| ) have been used and an identical proportionality factor is understood. The unperturbed contribution 
|£| 2 simply oscillates between values 0,4 with a period (tt/N). The function \ij\ 2 , instead, is strongly oscillating for 
increasing N, developing 2(N — 1) local minima and a sharply-peaked absolute maximum at uoAt = ir. Constructive 
interference is highest at the maximum, leading to a value Mj^x = 16 N 2 that can be very large, while destructive 
interference strongly damps the function for usAt < tt/2, which is a zero for both |77| 2 and |£| 2 for any N. One can 
show that 



\ri(N,ujAt)\ 2 < \£(N,ujAt)\ 2 on [0,tt/2] for any N . (44) 

Back to decoherence properties, condition (|44|) means that, for a mode at frequency u>, a finite region u>At < tt/2 
exists, where the contribution to decoherence is smaller in the presence of pulses. Since the "correcting region" is 
entered for small At values, this effect takes place in the regime of rapid flipping we expected. Since, moreover, 
smaller At values require longer pulse sequences in order to evolve the system over the same interval, it is useful to 
consider an interesting limiting case. 
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C. The limit of continuous flipping and suppression of decoherence 



Let us study an idealized situation represented by the following mathematical limit: 

At -> , 

N -> oo , (45) 
2NAt = t N -t ■ 

It is convenient to rewrite the decoherence function ( pl| ) by exploiting (|3|) to separate formally the unperturbed and 
the interference contributions. Thus 
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T P (N, At) = £ '^p 1 coth (jgL j 1 - /fc (iV, At) , (46) 
where 

At) = 2 6(At A } - V e2i(n-i)<-*At (47) 

In this way ro(tjv — to) is recovered by putting = for each mode, see (|l3|). We evaluate the asymptotic limit of 
fk(N,At) as follows: 

lim f k (N,At) = —r t — lim ^ '- V 2Ate Mfct - 1 

At-o v ; 1 - e l "<=(*«-*°) At->o At ^ 

n— 1 

lim i / (W^ s 



1 _ e iwfe(*jv-io) At— At 



to 



lim 

At— o 



sinwfeAt 1 — coscjfeAt 
+ « 



1 . (48) 



ujkAt wfeAt 

If this result holds for an arbitrary field mode, then the implications for the decoherence properties are transparent: 

lim o T P (N, At) = , (49) 

i.e. in the limit of continuous flipping, decoherence is completely and exactly eliminated for any temperature and any 
spectral density function. 

Obviously there is no hope that a continuous limit of this kind would be ever attained in practice. However, 
situation ( [45| ) should be approached if At is made small compared to the fastest characteristic time present in the 
dynamics of the system. From the considerations of Sec. II, the environment correlation time r c certainly provides 
a lower bound since there is no spectral content of the environmental noise at frequency higher than u) c . Hence, we 
expect that a sufficient condition in order to meet ( j45| ) is 

w c At<l. (50) 

The question on whether wc can do better than this may be not completely obvious, since time scales different from 
t c are also involved in the decoherence process. In what follows, we try to understand this point both by presenting 
a physical explanation for the decoherence suppression and by analyzing some specific situations. 

D. Physical interpretation 

People familiar with quantum Zeno effect f20|| may have found some similarities with the behavior we are discussing. 
In particular, the basic mathematical modification associated with pulses is a function \r)k\ 2 , which is 0(w|At 4 ) for 
short time intervals, compared to |£fc| 2 = 0(cj^Ai 2 ). Moreover, ( |45|) is formally the same continuous limit involved in 
many quantum Zeno related proposals, notably the one by Cook p3[ . In both cases, a preexisting dynamics - the qubit 
+ bath evolution here, a stimulated two-level transition in Cook's scheme - is modified through a pulsing procedure. 
Pulses respectively represent spin-flip interactions, short enough that the action of the bath is made negligible, or 
measurement pulses, long enough that the coupling to the external environment (the measuring apparatus) is made 
appreciable. A dynamical inhibition phenomenon occurs when pulses become sufficiently frequent. One can say that 
two opposite configurations are realized for decoherence: with continuous measurements, the interaction with the bath 
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is always "on" and internal dynamics becomes frozen; in the limit of continuous nipping, it is the two-level controlled 



dynamics that dominates, and the interaction with the bath tends to be always "off", as indicated by (49), However, 
the analogy stops from a more physical point of view. 

A more interesting interpretation of the decoherence suppression (f49| ) can be obtained by connecting it to effects 
already observed in magnetic resonance experiments, like spin-echoes, solid-echoes, or spin-flip narrowing All 
of these are basically time-reversal experiments. In order to capture the basic physical mechanism of our model, let 
us go back to the elementaty spin-cycle (Sec. IIIA) and look more carefully at the evolution operator (p3|). This is 
made of two pieces: a free evolution during the first At, followed by an evolution governed by ( |28|) during the second 
interval At. It is this transformed operator that simulates the effect of a time-reversal. In fact, backward propagation 
during the second part of the cycle would correspond to 



U tot (t + At, to) = exp [ - ^ (bW^'U&t) - b k e-^e k (At)) ] , 

^ k ' 



(51) 



By comparison with (|2^), forward propagation in the presence of rf- kicks only differs for an additional term e lWfcA * 
affecting each mode. This phase factor, due to the dynamics of the bath oscillators, is ultimately responsible for 
the decoherence properties in the pulsed evolution of the system. If not for this phase difference between (]2|) and 
(^T|), we would have 77^ (At) = and no decoherence. Instead, reversal is approximate since the bath restarts at time 
to + At after the first pulse with a dephased initial condition (see also Appendix). However, if e lWfcAt « 1 for each 



mode, then the couple of kicks produces an exact time-reversal and, by iteration, we arrive at (49). Equivalently, if 
the bath Hamiltonian can be considered as a constant, then the total Hamiltonian (|l|) acquires a minus sign and the 
system retraces the previous evolution. The validity of this condition depends on the time scale we are considering. 
For a single mode of frequency u>, the time needed to produce appreciable dephasing is t w so we expect 

decoherence correction for t/ At > 1. This is in agreement with both the interpretation of (|4|) and the semiclassical 
NMR argument pjj ], where motional effects are predicted for r comparable to the mean time spent in a given spin 
configuration (here At). For the whole environment, the correlation time r c is the minimum time scale over which the 
dynamics is approximately unchanged, and the same reasoning leads to a physical explanation of the rapid flipping 
condition (pT)|). 

As a final remark, we point out that the mechanism accomplishing time-reversal in our model is purely macroscopic 
in the sense that no reference is made to the dynamical state of the system. This is different from the familiar case 
of the Maxwell demon, that effectively reverses a dynamical evolution by operating over some microscopic variables 
(like velocities) at a given instant. Rather, the reversal is obtained by changing the sign of the system Hamiltonian 
through the action of suitable external fields (control). A different kind of demon, the so-called Loschmidt demon, 
has been introduced by some authors to account for this behavior p4| . In this terminology, the spin-flip procedure 
realizes in principle a Loschmidt demon for a decohering qubit. 



IV. ANALYSIS AND EXAMPLES 



In this section we try to give some semiquantitative picture of the decoherence mechanisms discussed so far. We 
focus our attention on a representative class of reservoirs, corresponding to so-called Ohmic environments. The 
appropriate spectral density is given by (|l6|) with n = 1. The time dependence of the decoherence is summarized by 
the following expression: 



r(t) 



Ct 



dto e 



1 



1 — cos ujt 



AT, At) 



(52) 



where we assume to = 0. Eq. ( |52| ) reproduces the unperturbed behavior of Sec. II when / = 0, in which case 
r(t) = ro(t), Eq. (fl3|). In the presence of N spin-flip cycles, Eq. j5^ ) is found from the continuous limit of Tp(N, At) 
in ([ll]), with f(uj 7 N 7 At) given by (|47|). According to Sect. Ill, in this case we are interested at decoherence after a 
complete pulse sequence, i.e. t = tjy = 2 A At. For a fixed strength a of the system- reservoir coupling, the properties 



of the environment enter (|52|) with two parameters, ui c 



T 



'0 



Let us first analyze decoherence in the absence of any correction, / = 0. Qualitatively different behaviors arise 
depending on the relationship between the cut-off frequency lo c and the thermal frequency cop = T . Typically, two 
extreme situations are considered. 

(i) lu c <C ujp: high-temperature limit or classical environment: 

Decoherence dynamics is relatively easy to describe in this regime since, due to the exponential dependence on the 
cut-off, uj c is actually the only characteristic frequency accessible to the system. The environment looks classical 
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in the sense that its quantized structure cannot be appreciated compared to the thermal quantum up. Accordingly, 
thermal fluctuations always dominate over vacuum ones and, after a short transient 0(rp) where decoherence is almost 
ineffective, dynamics becomes very fast and coherence is lost completely after a time comparable to r c . In general, 
one expects that in this limit equivalent results are obtained by assuming a heat bath of classical harmonic oscillators. 
A mixed quantum-classical derivation for a two-state open system is given for instance in p5| . 
(ii) u c 3> up: low-temperature limit or quantum environment: 

In this limit a more complex interplay between thermal and vacuum effects arises. Thermal fluctuations are only 
effective for t > Tp and, due to the exponential suppression of n(u,T), they are almost totally contributed by low 
frequency modes u ^ up. The effects of vacuum fluctuations dominate on an intermediate region r c < t < Tp, 
a nonvanishing contribution remaining, however, at longer time scales t > Tp. The frequency composition of the 
fluctuations is now less clear: while modes below the thermal threshold are still responsible for the long-time dynamics 
t > Tp, frequencies in the range up to u c mostly contribute for t < Tp but are also present beyond Tp. As a consequence, 
even at thermal time scales t > Tp 3> r c , high frequency modes contribute appreciably to the decoherence process and 
characteristic times 0(r c ) are still relevant in the underlying dynamics. 

Typical decoherence curves e _r ^ for the Ohmic environment are found by numerical integration of and 
are shown in Fig. 1 for two choices corresponding to high- and low-temperature limit, u c /T = 10~ 2 and u c /T = 
10 2 respectively. The partial contributions due to thermal and vacuum fluctuations are indicated separately where 
possible. In the low-temperature case, a quiet (t < r c ), a quantum (r c < t < Tp) and a thermal (t > Tp) regime are 
easily identified in the process, as indicated above and discussed in more detail by many authors [ppf29(1. In both 
configurations, the qubit coherence decays exponentially fast once the thermal regime is well established, 

e -r« M e -t/t th ; (53) 

for a suitable time constant tth- In models where the decoherence rate is constant in time, Eq. ( p3| ) is usually assumed 
as the definition of a typical decoherence time, tdec = Uh- In our case, since the whole behavior of T(t) is required for 
a complete knowledge of the decoherence dynamics, the definition of a characteristic time for loss of unitarity is less 
clear. Oversimplifying things, the situation can be summarized as follows: for both the high- and the low-temperature 
regimes, a characteristic time exists, indicating the departure of coherence from unity. This time is determined by 
the shortest between the two time scales T c ,Tp. Once this transient is over, the duration of the actual decoherence 
process is at least comparable to r c in the classical environment, and to Tp in the quantum one. If this information 
is used as an estimate for a characteristic decoherence time, we are lead to tdec ~ 0(r c /a) for the high-temperature 
limit and tdec ~ 0(rp/a) for the low-temperature limit respectively. For identical values of a and u c , decoherence 
occurs extremely faster in the former case, as expected on intuitive grounds. 

We now come back to examine how decoherence is improved in the presence of spin-flip cycles. With the interference 
contribution / restored in (|52|), we have calculated numerically the decoherence values obtained when a fixed time 
interval t is divided in an increasing number of cycles of duration 2At, i.e. 

A *=^V> 7V = l,2,...,A max . (54) 

The behavior of e~ r ^ as a function of the pulse frequency 1/At has been studied, and the procedure repeated for 
different representative times. The results for the high- and low-temperature reservoirs considered above are shown in 
Figs. 2, 3 respectively. The unperturbed values of decoherence at the appropriate times can be read from Fig. 1. We 
see that, as predicted by (|4|), decoherence-correction starts as soon as the region T c /At > 1 is entered. For times short 
enough that u c t < 2N, this may by even accomplished with a single cycle. In general, once this condition is fulfilled, 
no further reduction of At is demanded to evolve the system in a decoherence-free way. However, a warning also 
emerges from Fig. 3: if flipping is not frequent enough, not only does the correction effect disappear, but decoherence 
can be actually made worse compared to the one in the absence of pulses. The explanation of this behavior is rooted 
into the interference mechanism that builds up the correction factor (|30|): in a sense, decoherence can be subtracted 
almost completely from the frequency range uAt < 1 only at the expense of enhancing the decoherence contributions 
from modes outside that region. This intrinsic feature of the model is also relevant to understand why, despite of 
the differences existing between the high- and low-temperature decoherence properties, the same condition u c At < 1 
is required to prevent both thermal and vacuum noises. In fact, this condition comes quite natural for a classical 
environment, but one might at first wonder why a weaker condition upAt < 1 would not suffice for the quantum 
case, even if the frequencies u < top contain the fraction that mostly contributes at times of the order of tdec- The 
reason for this failure is the presence of vacuum fluctuations. By satisfying condition up At < 1, we do get rid of 
thermal dephasing, but we do not correct completely vacuum noise until u c At < 1. Precisely, we are missing modes 
of intermediate frequency up < u < u c that, although of minor importance at long times, may introduce amplified 
decoherence contributions if not properly corrected. 
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We conclude with a few comments on the relevance of our procedure for quantum information processing. While 
the idealized character of the model prevents us from a quantitative discussion of implementation criteria, we can 
compare with the principles underlying current quantum error-correction proposals. Essentially, these are schemes to 
encode redundantly information in such a way that it can be restored also when errors due to external sources have 
occurred. The syndrome-identification and the error-correction stages may be regarded as a feedback configuration: 
suitable measurement protocols are required both in conventional schemes [^|,^|-[l5| and in alternate techniques based 
on the quantum Zeno effect [l^JT7[ | , while conditional logic is exploited in the coherence-preserving routines proposed 
in fLSfl . In any case, error-correction methods have the effect of reducing the error rate per unit time and, in order 
to be effective, they must be repeated at time intervals At shorter than the typical decoherence time of the system, 
i.e. t dec I 'At > 1. In comparison to quantum error correction schemes, our procedure exhibits two fundamental 
differences: no ancillary bits are required to store the information; no measurements are performed. In principle, 
these could be advantageous features, since encoding would be more efficient and, by avoiding measurements, no slow- 
down of the computational speed would be introduced. In addition, rather than reducing the error rate, this method 
would suppress completely the error source provided the appropriate condition t c / At > 1 is fulfilled. From a more 
practical perspective, it is the accessibility of this rapid flipping limit, demanding fast and short pulses, t c / At > 1 and 
Tp/At <C 1, that determines the viability of the procedure itself. If such req uirements can be satisfied, our method 
might be valuable in configurations where tdec tends to be shorter than r c [p6| or, even in case td ec is longer compared 
to r c , for systems where tipping the state is easier than exploiting conventional error-correction protocols. While the 
existence of an interaction able to implement a NOT gate by inverting the state is demanded for any two-level system 
relevant to quantum computation, both the present technological capabilities for applying 7r-pulses and the relevant 
environmental cut-off frequencies depend considerably on the specific physical system and the mechanism responsible 
for decoherence. Some important time scales for various prospective qubits can be found in M . 



Our work demonstrates the possibility to modify the evolution of a quantum open system through the application 
of an external controllable interaction. A prototype situation involving a two-level system coupled to a quantized 
reservoir in thermal equilibrium has been worked out in detail and dynamical suppression of quantum decoherence has 
been evidenced. From the perspective of quantum information, the analysis suggests a different direction compared 
to conventional quantum error-correction techniques, based on the idea of forcing the system into a dynamics that 
disturbs the decoherence process. Our present study for a specific example brings up, among other issues, the question 
of whether similar decoherence correction mechanisms would be operating under more general conditions, including 
either different open system dynamics, or different control configurations, or both. In particular, an interesting 
possibility could emerge from examining decoherence properties within a fully quantum mechanical description where 
the control degrees of freedom are explicitly included and the system is driven by a quantum controller as recently 
proposed in JjJ. 
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Compared to the interaction picture, the Heisenberg representation has two advantages: first, it does not require 
preliminary transformations on the state vector; second, it gives to a certain extent a more intuitive description of 
the spin motion. In this Appendix, we outline the evolution of the qubit coherence, by restricting ourselves to the 
first elementary spin-flip cycle. In the Heisenberg picture, the relevant information is contained in 



since, by averaging over the quantum state, (a+(t)) = pm{t). As in Sec. Ill, we evaluate the qubit dynamics under 
the separate action of the spin-bath Hamiltonian Hsb and the radiofrequency perturbation H r f. The description of 



V. CONCLUSIONS 
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APPENDIX A: HEISENBERG REPRESENTATION 




(Al) 



11 



a 7r-pulse turns out to be extremely simple in the Heisenberg representation. Nothing happens to the bath operators 
b k ,b k in the limit of instantaneous pulses, while spin dynamics is governed by the equations 



o a , (h s + H rf (uj ,t) 



a = x,y,z 



(A2) 



with H r f(u>o,t) given in ([L9|). By denoting with t P ^ the instants immediately before (after) a pulse respectively, a 
very simple result is found: 



<?z(tp) = - <Jz{t P ) , 
a+(t+) = (o-+(tp))t . 



(A3) 



The action on a z corresponds, in particular, to the pictorial spin-flip effect operated by a 7r-pulse. Now denote as 
Gtot{U,tj) the operator evolving coherence from ti to tj in the absence of rf-pulses, i.e. 



a +(tj) = Gtot(U,tj) a+(ti) . 



(A4) 



Then, using relation (AS) for a+ twice, we find the following representation for the coherence evolution during the 
first complete cycle: 



to be compared with 



a+(t a + 2 At) = G tot (t , t + At) a+(t ) G\ ot {t a + At, t + 2 At) 



a+ (t Q + 2 At) = dot (to + At,t Q + 2 At) G tot (t , t Q + At) a+ (t 



(A5) 



(A6) 



in the absence of pulses. Even before knowing the explicit form of Gtot, we see from ( A5)-([Aq) that the presence of 
a time-reversed evolution during the second interval At of the cycle is already enucleated at this stage. 

In order to evaluate the propagator Gtot, the Heisenberg equations for the coupled spin + bath motion have to be 
solved. From Hamiltonian (Tj]) we get 



( b k = -iuj k b k - ig k a z , 

b\ = +iuj k b{ + ig* k a z , 

cr + = iu (T + + 2iY, k (9kb\ +g* k b k ) <r + 

o z = 0. 



(A7) 



Since instantaneous pulses introduce discontinuous changes in operators, the propagators Gtot(to,to + Af), Gtot(to + 

At, to + 2 At) have to be considered separately, by solving (A7) with initial conditions at t = to, t = t~p = to + At 

respectively. 

(i) to^t + At: 

Since (J z (t) = &z(to), the equations for the bath variables are completely solved by 



b k (t) = e-^ t - t ^b k (to)-aSo)—(l-e- i ^ t - t '>A , 

u) k \ ) 

and b\(t) = (b k (t)y . These solutions should be inserted in the expression for cr+(t): 



(A8) 



<j + (t) = Texp <^ i / ds ujq + 2(g k b\(s) + g* k b k (s)) \ <J+(t ) 



(A9) 



The time-ordered exponential can be evaluated exactly and the following result is found for the first propagator: 



Gtot (*o, to + At) = exp <^ iuj At + At - 



,.-\9k 



sin ui k At 



(l-a z (t )) 



CX P { - E ( & !(*o)6(At) - b k (t Q )C k (At)) } , 



(A10) 



where the same notation for £ k (At) has been used, Eq. (]?]). 
(ii) to + At ^ t + 2At: 
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By exploiting (A10), we can immediately write down the expression for the propagator G\ ot (to + At, to + 2At) in 

sinw^Ai v 



terms of the new initial condition at t = tp. 



Gl ot {t + At, t + 2At) = exp <j - iuj At - U l -^(At - — Vl - a z (t+)) 

LO k \ Ul k t 



exp \ + (b{(t+)^(At) - b k (t+)C k (At) 



(All) 



Now everything can be evaluated with respect to the initial time of the cycle. We exploit (A3) for <r z (tp) and, since 
bk{tp) — bk(tp) — b k (t + At), bath operators evolve as 



b k (t + At) = e-^'bkito) - tr z (to)— ( 1 - e' 



(A12) 



and Hermitian conjugate. As we discussed in the text, the presence in (A12) of an initial condition dephased by 
e ±iuj k At £ or eacn environmental mode is the ultimate source of decoherence. We arrive at the following expression for 
the cycle evolution: 

<7+(t + 2At) = expiitpi(At)(l - <7 z (t )) + ^ (&l(*o)%( A *) - b k (t )v* k { At )) + m( At ) 



a+(t )exp< ia z (t )<P2(At) - ipi(At)(l + a z (t )) > , 



where 



(At) = 4 V^fA< 

^ LU k \ 



sinwfeAt 



L0 k 



la I 2 

f 2 {At) = 8^-^-sinw fc Ai(l - cos uj k At) 



k ^ 



The final step is to calculate the coherence evolution as 

poi(*o + 2At) - ((T + (t + 2At)> = 01 Tl 's{ (Pb ® Ps) <r+(t a + 2At) } \j) . 

3=0,1 



By inserting (A13)-(A14), phase factors drop out and we find 

Poi(io + 2Ai) = p i(io)c- r - (w=1 ' At) 



(A13) 



(A14) 



(A15) 



(A16) 



in agreement with the result found from Eqs. (f4(i|)-(|4l|) in the Schrodinger representation. The procedure can be 
generalized to arbitrary N and the complete expression Tp(N, At) is thereby recovered. 
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FIG. 1. Qubit decoherence as a function of time for an Ohmic environment, Eq. (52) (/ = 0). Time is in units of T _1 
and the values lo c = 100, a = 0.25 have been chosen. High- and low-temperature behaviors are shown, (H) uj c /T — 10~ 2 and 
(L) uj c /T — 10 2 respectively. The contributions arising from the separate integration of thermal and vacuum fluctuations are 
displayed in the latter case, e~ r(t) = e~ r * h(t) • e~ r " (t) . 

FIG. 2. Qubit decoherence in the presence of rf-pulses for the high-temperature configuration, uj c /T = 10~ 2 . For a fixed 
time, each point corresponds to a number N of cycles, N = 1, . . . , -/V ma x = 10. The results from Eq. (52) are plotted as a 
function of the normalized pulse frequency T c /At. The unperturbed values of decoherence are read from Fig. 1 (H). 

FIG. 3. Same as in Fig. 2 for the low-temperature configuration, lo c /T — 10 2 . The maximum number of spin-cycles is equal 
to iV max = 30 in the simulations at uj c t — 1.0, 10, while -/V max = 100 at uj c t = 10 2 . The unperturbed values of decoherence are 
read from Fig. 1 (L). 
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